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Abstract 

We study the graviton self-energy function in a general gauge, using a hard 

thermal loop expansion which includes terms proportional to T 4 , T 2 and 

log(T). We verify explicitly the gauge independence of the leading T 4 term 

and obtain a compact expression for the sub-leading T 2 contribution. It is 

shown that the logarithmic term has the same structure as the ultraviolet pole 

part of the T = self-energy function. We argue that the gauge-dependent 

part of the T 2 contribution is effectively canceled in the dispersion relations 

of the graviton plasma, and present the solutions of these equations. 
ll.10.Wx 
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I. INTRODUCTION 



When the temperature T is high compared with the typical momentum scale but well 
below the Planck scale, all the n-graviton thermal Green functions can be computed in the 
one-loop approximation using the hard thermal loop expansion. There have been many 
investigations where this approach has been employed [0-0]. An important property which 
is now well established is the gauge invariance of the leading high temperature contributions 
of all n-graviton thermal Green functions. The explicit from of these contributions can be 
obtained using the equivalence which exists between the formalism of Boltzmann transport 
equation and the high temperature limit of the thermal Green functions in field theory |§. 
Using this approach (which is explicitly gauge invariant) one can easily show that the leading 
part of all n-point one-loop thermal Green functions is proportional to T 4 [|§ . These results 
have also been obtained by standard Feynman diagrammatic calculation in the Feynman- 
DeDonder gauge for the one- and two-graviton functions H as well as for the three-graviton 
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One of the interesting physical applications of the one- and two-graviton functions is the 
study of the dispersion relations |4j which follow from linear response theory . Since the 
relevant physical quantities are obtained from the poles of the propagator, it is important 
to verify the gauge independence of this procedure. While this is automatically fulfilled 
by the leading high temperature contributions, the inherently gauge dependent sub-leading 
contributions to the thermal Green functions require a more detailed investigation. A 
similar situation occurs in the case the Yang-Mills theory where it is known that a gauge 
independent set of dispersion relations can be obtained from the gauge dependent thermal 
two-gluon function 0,0. As far as we know, in contrast with the case of the Yang-Mills 
theory, a gauge independence proof of the dispersion relations in quantum gravity beyond 
the leading order is still missing. 

The purpose of the present paper is to investigate this problem in the case of gravity us- 
ing the standard Feynman diagrammatic approach. We will compute the 1- and 2-graviton 
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functions up to sub-leading contributions, in a class of general gauges. (Some intrinsically 
gauge independent sub-leading contributions have been considered before using the theory 
of a scalar field in interaction with gravity ||.) We employ the imaginary time formalism 
|T3| and express the one-loop thermal Green functions in terms of on-shell forward scat- 
tering amplitudes (the "Barton amplitude") JTJJ], properly generalized in order to account 
for the quadratic denominators which arises in the free graviton propagator when a gen- 
eral gauge fixing term is employed |L5| (see also appendix [A]). This approach enable us to 
explore some of the general properties of the exact graviton self-energy without having to 
carry out explicitly the non-trivial spatial momentum integrations. It is also much more 
straightforward to perform the hard thermal loop expansion when we start from the forward 
scattering amplitudes. Using this approach we were able to obtain explicit results for the 
T 2 and logarithmic contributions to the graviton self-energy. 

Unlike the leading high temperature terms, for which the gauge independence is con- 
firmed by our calculation, the sub- leading contributions are gauge dependent. These con- 
tributions will be employed in the study of the dispersion relations for the transverse and 
traceless gravitational modes ||. 

This paper is organized as follows: In Sect. [TI| we present the Lagrangian and the basic 
definition of the graviton field from which the Feynman rules are derived. We also discuss the 



identities which follows from the gauge invariance of the theory. In Sect. |TTT] the main results 
of the calculation of the one- and two-graviton functions up to the logarithmic contributions 
is presented. A very compact expression for the T 2 contribution is obtained and we verify 
that the logarithmic contribution is proportional to the ultraviolet pole part of the T = two- 
graviton function. We also derive the general transformation of the two-graviton function 
under a change of graviton representation. In Sect. [IV] we verify explicitly that the gauge 
dependence of the sub-leading T 2 contributions to the two-graviton function is effectively 
canceled in the dispersion relations. We then present the solutions of these equations, which 
include corrections of order T 2 to the leading T 4 contributions, describing the physical modes 
for the propagation of waves in a graviton plasma. 



II. FEYNMAN RULES AND IDENTITIES 



The Feynman rules for the graviton propagator and self-interactions vertices are obtained 
from the following underlying Lagrangian 

9 1 r) I nnV LV 

c = -i^9 R+ -r^ (ppV^sT) (doyFgsT) + v 5 J v x ; * = Vs^G, (2.1) 

where R is the Ricci scalar, G is the Newton constant and the parameter £ defines a family 
of gauges. (£ = 1 is the Feynman gauge and £ = is the Landau gauge). The quantities Xu 
and f] x are the ghost fields and the function e(x) is the infinitesimal generator of coordinate 
(gauge) transformations 

x» -^x" + e^(x). (2.2) 

The calculations in quantum gravity are conveniently performed using the graviton field h^ u 
defined in terms of the tensor g^ u as 

yf^g^ = TT + Khr, (2.3) 

where r] fJ-l/ is the Minkowski metric. 

The Feynman rules can be obtained in a straightforward way substituting ( |2.3| ) into ( |2.1|) 
and performing a perturbative expansion in k. The Oth order terms are quadratic in the 
graviton field and yield the following expression for the graviton propagator 



(1 — 



(2.4) 



In the appendix [FJ we give all the other relevant Feynman rules employed in this work. 
The choice of the graviton field parametrization given by (|2.3|) restricts the gauge pa- 



rameter dependence only to the propagator (|2.4|) , since in this case the second term in ( pTT 



is exactly quadratic in the graviton field h. This is similar to the general linear gauges 
in Yang-Mills theories. Therefore, the gauge dependence of the Green functions computed 
from these Feynman rules can be traced back to Eq. ( |2.4| ). 

The leading high temperature contributions of all one-particle irreducible thermal Green 
functions are related to each other through tree-like Ward identities in the same way as the 
basic tree vertices |||5|||. These hard thermal loop identities have been verified for both 
Yang-Mills theories and gravity and generalized to any gauge theory whose generators form 
a closed algebra ||. For our present purposes it will be sufficient to consider the identity 
involving the two-point function. A simple example is provided by the following tree Ward 
identity, arising from the invariance of the pure Einstein action Sq = 2/k 2 / d A Xy/—gR under 
the transformation given by equation (|2.2| ), 

X^ UiX (k)S^^(k) = 0, (2.5) 

where 

X uJxi k ) = k ^VuX + - k x r]^, (2.6) 



is the tensor generated from the transformation of the graviton field under ( |2.2p and 

k 2 

•S^iii ^ 2i / 2 (^) = w ( r ?MiM2 ? 7^i^2 + Vm^Vvi^ ~ V^iviV^^) 

~\~ g { k fJ,l k fJ.2VuiU2 ~l~ k fJ,l k U2VvifJ,2 ^Vl ^2^1/^2 kuik^VfJ-l^) 



(2.7) 



comes from the quadratic term in the action without the gauge fixing term (it is the inverse 
of the propagator in the limit £ — ■> oo). 

The tree-like identity which holds for the high temperature limit of the two-graviton 
function would be identical to Eq. ( |2.5| ) if the one-graviton function (the tadpole), shown 
in the diagrams in Fig. 1, were zero. The modification introduced by the tadpole changes 
the right hand side of Eq. Q2.5p to a non zero quantity when Sj®^ M2 „ 2 (^) is replaced by the 
leading high temperature contribution of IP 1 1/1 M2!/2 (fc, u), given by the diagrams in Fig. 2 
(u is a time-like normalized four- vector representing the local rest frame of the plasma). 
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This contrasts with the analogous situation in the case of Yang-Mills theories where the 
anti-symmetry of the group structure constants trivially makes the tadpole to vanish. As a 
consequence of the non- vanishing tadpole, the general BRST identities will not hold for the 
exact finite temperature graviton self-energy However, as we will see in the next section, the 
tadpole diagrams can be computed exactly, yielding a result proportional to T 4 . Therefore, 
if we split ^ 2U2 (k, u) as 

it— = nSf(M) + n£? M2 "(M), (2.8) 

the BRST identities derived in the appendix |C| will hold for the sub-leading contributions 
IT^b 1 fJ,2U2 (k,u), so that the following identity is satisfied 

X^ x (k)UZ^" 21/2 (k,u)X^ 2S (k) = 0. (2.9) 



This identity is analogous to k^k u Il^ D = 0, where IT^f D is the exact gluon self-energy |12 



Since all the gauge parameter dependence is restricted to the sub-leading contributions, 
these identities have an important role in the cancellation of the gauge dependence in the 
dispersion relations. 



III. THE ONE- AND TWO-GRAVITON FUNCTIONS IN A GENERAL GAUGE 

In this section we will present the details of the calculation of the one- and two-graviton 
functions. Let us first consider the contributions from the two tadpole diagrams in figure 
1. The most involved diagram is the one shown in figure la, since both the 3-graviton 
vertex and the general gauge propagator are involved. Using Eq. ( |B4| ) and the propagator 
( |2.4| ), the straightforward contraction of indices yields a result which is independent of the 
parameter (1 — £). Therefore, the resulting expression is identical to what is obtained in 
the Feynman-DeDonder gauge involving only the usual quadratic denominators. Using the 
Eq. ( |XT| ) in the simple case when i = 1 and j = the following result for the one-graviton 
function is readily obtained 



6 



\q\ 


2 d\q 


1 


2 


1? 


e \q\Q-u/T _ 1 



(3.1) 

= -«p y = «| - 4u M u„) ; P = ^T A 

where / dQ denotes the angular integral and the four vector is on-shell with components 
given by = (l,q/\q\). 

The diagrams contributing to IP m 1121/2 (k, u) are shown in figure 2. The contributions 
associated with each of these diagrams will involve integrals like the one shown in Eq. ( |A6|) . 
From the structure of the graviton propagator given by Eq. ( |2.4| ) we can see that the diagram 
in figure 2a is such that each of its terms will involve integrals with i, j = 1,2, while in the 
case of the diagram shown in figure 2b, all the corresponding integrals have the form of 
the first term of Eq. ( |A6[ ) with j — and i — 1, 2. In the case of the ghost loop diagram 
shown in figure 2c, all the terms will involve integrals with i = j = 1. Let us first consider 
the leading high temperature behavior of these integrals. In this limit, we can perform a 
hard thermal loop expansion of the integrand such that the terms with i,j > 1 will all be 
sub-leading. For the terms with i,j = l we use expansions like 



11 lk 2 1 {k 2 ^ 2 
q2=Q ~ 2q~k~\J^f + 8j^. k ) 



1 + ^ + -' (3-2) 



(q + k) 

The case i — 1, j — (from the diagram in figure 2b) is similar to the tadpole diagram 
giving an exact T 4 contribution. In this way we obtain the following result for the leading 
behavior of the graviton self-energy 



Yfleading ft, „.\ _ JlP ( fMiQ viQ fi 2 Q ' V2 . QjH Q ^2 Qu 2 , 

^u m ^U) - K -J — I — + — + 



Qui Q fii Qu2 ^ ^uiQu\ Qn-iQ^x k QuiQviQ \iiQvi 



(3.3) 



k-Q k-Q (k-Q) 2 J 

It is worth mentioning that though a naive power counting would allow for a gauge parameter 

dependence from the third term in the second line of Eq. ( |2.4j ), the final result ( |3.3| ) is gauge 

independent as one would expect on more general grounds |3|. Combining the Eqs. 

fl2.9|) and (|3.3| ) we obtain the following identity for the exact self-energy 



(3.4) 



where in the last term we have used ( p.l|) . Since the integrand in the right hand side of 
the above expression is an elementary expression without denominators, the same should be 
true for its left hand side, up to terms which would vanish after integration. Our calculation 
shows that, in fact, the expression obtained from the diagrams in figure 2 is such that the 
exact integrand of X^ Vi A IP 1 1/1 ^ 2U2 X^ 2 \ 2 s does not involve any denominators, being identical 



to the integrand in the right hand side of ( |3.4| ). 

We have extended the hard thermal loop expansion in order to obtain the T 2 and the 
logarithmic contributions, which are yielded respectively by the terms of degree —1 and 
—3 in \q\ (terms of degree —2 in \q\ are absent due to the symmetry q «-> —q) from the 
expansion of the integrand in expressions like ( |A6| ) in the region of large values of \q\. After 
a long computation we have been able to find the following compact expression for the T 2 
contribution 

k 2 T 2 



n 



P\v\P2V2 



(k,u) 
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i Q(0) tP<tXS q(0) 
V1P2V2 \i\V\po ■ F °\Sp2V2 



(3.5) 



where 



jpaXS / v *" " 

i F = / — 



Arc 



AQ p Q a n xs + AQ x Q s r]P a - ZQvQ^ 5 k 2 QfQ a Q x Q 5 



(k-Qy 



(k-QY 



(3.6) 



and 



jpaXS 

L G = 



Air 



Q p Q a r] xs + Q x Q d n 



X/^jS^pa 



Using Eq. 



(k-Q) 2 

and the structure of (|3.5| ) we immediately conclude that 



(3.7) 



X^^\k)Ull^ V2 (k,u) = 0. 



(3.8) 



This result can be understood in the context of the BRST identities, using the results of 
Appendix y. It is remarkable that though this T 2 contribution is gauge dependent, it is 
transversal to X^^ 1 A . 
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It is straightforward to obtain the explicit results for the angular integrals in Eqs. Q3.3D , 
( [3.6| ) and ( |3.7| ) in terms of a tensor basis such as the one shown in the table |. Using the 
following decomposition 

r 14 

/ dQf^^i^Q) =Y,Ci(k,u)Tt lVlll2U2 , (3.9) 
•* i=i 

where fvi v u x * v *(k ) Q) is a function of degree 2 or in Q respectively for the leading T 4 or the 
T 2 contributions, the coefficients Ci(k, u) are obtained contracting both sides of Eq. ( |3.9| ) 
with each of the 14 tensors of table |. The solution of the resulting linear system of 14 
equations is given in terms of integrals like / dQ (k ■ Q) r , which can be easily evaluated. 

In the case of the logarithmic contributions the resulting angular integrals / dQ can all 
be parametrized in a Lorentz covariant way in terms of the 5 tensors T^ 1 " 1 ^ 2 " 2 ( j^i"^ 1 ^ 
jiMi^iM2f2^ rppivi/x2V2 ^ j^i^iM2f2_ rjv^ resu [^ can expressed in terms of the T = graviton 
self-energy |l6j in the following way 



njTf W2 (M) = iog(T)n^^(fc), (3.io) 

where U^ iUifl2U2 (k) is the residue of the ultraviolet divergent T = contribution which is 
obtained from the calculation in n = 4 — 2 e dimensions. The fact that both the log(T) and 
the ultraviolet divergent contributions have the same structure has been also verified for 
the two and four-point functions in QED |l7j and for the two- and three-point functions in 
Yang-Mills theories ||18|| . These results are special examples of the rather general arguments 



presented in |19| . 



From the results for the thermal one- and two graviton functions we can write the 
following expression for the thermal effective action 

Sterm[g}=r, u h^(0)+ J d*k h"™ (k)Up lvl „ 2 „ 2 (k)h^ (~k) + • • • . (3.11) 



Here we will use Eq. ( 3.11|) in order to derive the expressions for the new one- and two- 



graviton functions which arise when one uses the graviton representation 

g^ = V^ + ^9^u- (3.12) 



The corresponding expressions will be employed in the analysis performed in the next section. 
Using Eqs. ( |2.3| ) and ( |3.12| ) one obtains the following relation for the graviton fields in the 



two representations 



(3.13) 



Inserting Eq. ( |3.13|) into Eq. (|3.11|) and using the traceless property of T^ u [cf. Eq. (|3.1| )1, 
we obtain 

Stermb) = ^,^(0) + J d 4 k5 9 ^(k) ^ (k) S 9 ^(-k) +■■■, (3.14) 

where 

fv = -r^ (3.15) 

and 

n" 1 " 1 " 2 ^ (k,u) = (k,u) ~ (u^ ia a 7]^ U2 + w^ a a rf liVi - ^ir/zj^V 21 ' 2 

-(-PWlflW 1 ^ _|_ pM2"2jyW I 'l 
_PM1A«2^1'1^2 _ p^l^2^MlM2 
pi'iA'2^A 1 i^2 YL llU2 rj l ' lfJ ' 2 ^J 

(3.16) 



We remark that while the derivation of Eq. ( |3.16| ) is rather simple and general, a direct 



calculation of nw" lM21 ' 2 (k), on the other hand, would involve the manipulation of more 
complicated Feynman rules where the gauge fixing term from Eq. (|2.1|) would contribute to 
all the n-graviton vertices. 

IV. THE GRAVITON DISPERSION RELATIONS 

The thermal graviton self-energy has been employed in order to investigate the propaga- 
tion of gravitational waves in a plasma ||. This can be done studying the poles of the full 
graviton propagator (dispersion relations) which is obtained from the effective action 
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S[g] = S G [g] + S flx [g} + S term [g} (4.1) 

where Sc[g] is the Einstein action, Sf ix [g] is the gauge fixing term and S term [g] is given by 
Eq. ( |3.11| ). In this section we shall apply the results for the graviton self-energy up to the 
sub-leading T 2 contributions in order to investigate the gauge dependence of the dispersion 
relations. 

Since the tadpole contribution to St erm \g\ yield a non-zero energy- momentum tensor in 
the Einstein equation 

^1 = 0, (4.2) 

a self-consistent calculation of the full graviton propagator has to take into account a curved 
background so that 

9w = ff£) + &9v» ( 4 -3) 

where gftj is the solution of the Einstein equation (|4.2|) and 5g^ v is the metric fluctuation. 
From the corresponding second order variation of the effective action 

S 2 S[g] = -\j d*x^gW Sg^P^^Sg^, (4.4) 

the graviton propagator can be obtained taking the inverse of p^ Ul i 12U2 , 

The contributions to p^ lUl ^ 2U2 from the first two terms in (|4.1| ) are well known |I, |20 , 21 



They involve components of the Riemann and Ricci tensors and the scalar curvature. Re- 
stricting the analysis to a metric background which is conformally flat, the components of 
Riemann tensor can be expressed only in terms of the Ricci tensor and the scalar curvature. 
Since Eq. Q3.1|) yields a traceless energy-momentum tensor, the Einstein equation ( |4.2| ) 
(with vanishing cosmological constant) implies that the scalar curvature is zero and that 
the Ricci tensor is proportional to Eq. ( |3.1| ). Using geodesic normal coordinates the thermal 
contributions to p^ 1 " 1 ^ 2l/2 can be obtained from Eq. ( |3.14| ) After a straightforward tensor 
algebra one obtains the following expression 
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+- ( 1 - — ] rf^rf^k 2 - -ri^^ri^k 2 
2 \ 24/ 2 

-8tt Gp (g^ 1 " 2 (4u w ^ 2 - r/^ 2 ) + ^ lUl (Au» 2 u V2 - r)" 2U2 )) > ( 4 - 5 ) 
+ (symmetrization under /i 2 <-> ^2) 

+ (symmetrization under /ix <-> z^) j — 167rG'n Ml ^ 1 ' 12!y2 (k, u) 

where IP 1 " 1 1X21,2 is given by Eq. ( |3.16| ) with the leading and the sub-leading high-temperature 
contributions from rP liyiAt2 ^ 2 given respectively by Eqs. (|3.3|) and fl3.5|) . 

Because of the coordinate invariance of the problem we have to impose physical con- 
straints on the metric fluctuations. The imposition that the spin one and spin zero degrees 
of freedom do not propagate constraints the metric perturbations Sg^ to be transverse 
and traceless, respectively ||. These conditions imply that we only have to consider the 
transverse and traceless components of (p a ^^ u \ in the linear response equation 

§g a/3 = -lQnG (p a ^A _1 ST"". (4.6) 



An explicit basis of TT-tensors can be found imposing the TT conditions on a general linear 
combination such as the one on the right hand side of Eq. ( |3.9| ). It is also convenient to 
choose these tensors as being idempotent and orthogonal to each other. This leads to the 
following set of TT-tensors 
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Tf lVl ^ V2 = ciT^ u ' » 2U2 , I = A,B, C. (4.7) 

i=l 

where the coefficients en are given in the table 0. 

This result is in agreement with the one obtained by Rebhan in reference 0] (except for 
a small misprint in the 9th line of the 2nd row in table 2). As a simple checkup of this result 
we note that at zero temperature there is only one TT-tensor given by 

q-fllUx fl 2 ^2 _ q-^lVl M2^2 _|_ ^Ml^l M2^2 _|_ ^"Ml^l M2^2 ^ 

so that in the table [TT]the lines 1, 4, 8, 12 and 13 (which gives the coefficients of the Lorentz 
covariant tensors in the table []]) must add to a Lorentz scalar (or a pure number) and all 
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the other lines must add to zero. Once we know a certain set of coefficients q in the basis 
given by table | and the explicit form of the TT-tensors given by Eq. (|4.7|) a straightforward 
calculation gives the following result for the coefficients in the basis of TT-tensors 

c A = 2ci 



/ k 2 -(k-uf \ 
C ^ 2 ( Cl + ^"^ C2 j • (4.9) 



cc 



( Ak 2 -(k-uf l(k 2 -(k-u) 2 ) 2 
* + 3 k 2 02 + 3 Wf C3 



It is interesting to note that the sub-leading contributions to the graviton self-energy are 
independent of the graviton representation. This property is satisfied because these con- 
tributions to anc [ Tpi^iM2^2 h ave the same TT components. Indeed, we see from 
Eq. (|3.16|) that apart from the tadpole contributions which have an exact T 4 behavior. The 



terms involving traces of IP 1 " 1 '* 21 ' 2 are either proportional to rj^ 1 ^ or rj^ 2U2 or both. Such 
terms have no components along any of the first 3 tensors of table |, so that they give no 
contribution to any of the coefficients Cj, i — 1, 2, 3 which appear in Eqs. ( f4.9Q . 

We have now all the basic quantities which are needed in order to express p^ 1 ^ 1 ^ 2 in 
the basis of TT-tensors as follows 

14 

pu-wrw* = CA T^ Ul ^ 2 + c B T£ lUlfX2l/2 + c c T^ Vl ^ 2 + Y / c l Tt lUl " 2U2 . (4.10) 

i=4 

The inverse 

14 

(P^ 1 1 ^ 2 2 ) = &aTa ii] vt, fj, 2 V2 ~l~ ^bTa mi v-j &cTa mi uj ti?,v?, ~t~ djTj M1 Vl (4.11) 

i=4 

can be determined from the relation 

(p^v lP *yl ppaw* = I + S^ 2 ) . (4.12) 

Using the transversality and idempotence of 7^"°^ as well as the identities 
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^A P a T i aaP = (z = 4---14) 



1 B pa l i 



for i ^ 6 
7^ for i = 6 



(4.13) 



rypV rppCT 0.(5 r\ 

J C pa I i — U 



4,8, 10- • -14) 



7f „ Tf aap = (i = 4 • • • 14; I = A, B, C) 



we obtain the following result 



1 



d A = — 

ca 



di 



1 

— ( 

2 



I " A ^ r r^ u rjipaafirr 
1 — ^"6 Cq i 6 pcr J 6 J-Bp.va.0 



(4.14) 



dc = — 1 



E.J „ rp^V rppoafiq- 
a i c j- L i pa 1 ] 2 Cpva/3 



(ij=5,6,7,9) 

We can now investigate the poles of the TT components of the propagator from the 
solution of the equations cj = 0, I = A,B,C. Using the Eqs. (fOj) with ci, C2 and C3 
determined from the decomposition of (|4.5f) in the basis of table |, the equations associated 
with the modes A, B and C can be written respectively in the form 
16tt G> 



k 2 



l-^Gf(l-()L 
15 V w 



5 1 4r 1 "A 8 A; 2 / 2 _ 1 4r 1 2 A 

- + - r A L r 2 + — — \8r L + -r A L r 2 - 4 

9 2 6 J tt 2 T 2 V 2 6 



16% Gp 



647T 

l + ^GT 2 (l-0 
15 



2 o4r 2 2 10 1\ 8A; 2 / r2r or4 2 2 

- - 2 r A L + - r 2 H + -=-=■ -5 r 2 L —2Lr +-r 

9 3 9 r 2 tt 2 T 2 V 3 



1-^Gf(l-()L 



8 , „ 28 1 \ 8 k 2 / 9 



r 4 _ r 2 



15 



(4.15) 
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where 



2 _ " , _ ko , k + \k\ , 

r = ¥~k ] L = m\ g k^m~ l (416) 



The Eqs. ( |4.15| ) reduces to the Eqs. (6.2) of reference in the special case when the 
sub-leading terms proportional to G T 2 or k 2 /T 2 are neglected. 



In view of the constraints imposed by the important condition (|3.8| ), the gauge dependent 
denominators in Eqs. (|4.15|) have a very simple structure. Since we assume that GT 2 C 1, 
the momentum-independent denominators can be expanded perturbatively. We thus see 
that all the gauge dependent T 2 sub-leading contributions to the dispersion relations give 
effectively corrections of order G 2 , which are of the same order as the two-loop contributions 
which we have disregarded. Hence, we conclude that to one-loop order, the dispersion 
relations are effectively gauge independent. 

The solution of the one-loop dispersion relations can therefore be obtained from Eqs. 
( [4.15| ), by setting the denominators equal to one. These solutions have been obtained in 
the reference Q in the leading high temperature approximation. In order to illustrate the 
magnitude of the T 2 sub-leading contributions let us consider the solution of Eqs. Q4.15| ) 
for real values of k$ and k. This corresponds to the propagation of waves supported by 
the graviton plasma. In the figure 3, where Q = ko/ (IQirGp) 1 ^ 2 and k = \k\/ (WirGp) 1 ^ 2 , 
we show the numerical solutions corresponding to the modes A, B and C and compare the 
leading results with the contributions which include the sub-leading T 2 corrections. We can 
see from these diagrams that for all TT-modes, the dispersion curves begin at a common 
plasma frequency cu p i and become asymptotically parallel to the light cone. 

The behavior of the dispersion relation can be determined analytically in the limiting 
cases of very small and very large momenta. When k — > 0, the common form of the dispersion 
relations is given by 

-§A l ~ !r Gr2 ) < 4 - 17 > 

For large momenta such that k 3> 1, the asymptotic forms of the dispersion relations become 
respectively 
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uj A = k A + 



UJ B = K B 



u c = k c 




) 



(4.18) 



The small wavelength limit given by the above relations can be understood by noticing that 
in this case one probes the plasma at small distances, where the medium effects on the free 
dispersion relations are relatively unimportant. On the other hand, for long wavelengths Eq. 
( |4.171 ) gives a substantial modification of the free dispersion relation due to the collective 
phenomena in the plasma. However, this modification is sensitive to the curvature of space 
which we have neglected in our analysis. (Such an effect yields corrections of magnitude 
(GT 4 )(GT 2 ), which are formally of the same order as the two-loop contributions). This is 
an interesting issue which deserves further study. 
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APPENDIX A: 



Here we show explicitly how to extend the method of Barton amplitudes in order to 
account for the contributions which arises from the quadratic denominators in the general 
gauge free propagator. We illustrate this technique by considering the following integral 



d 3 g 



1+6 dz 
s 2ni 



N(z) 



t(q;p) 



- + (z 



<-> —z) 



(Al) 



(q-q) (p-p) J 

where p = q + k, k = 2m7ri, m = 0, ±1, ±2 ■ ■ ■, q = (z, q) and i, j = 1, 2. This is the most 
general kind of integral which contributes to the two-point function when the imaginary 
time formalism is employed. The generalization to higher point functions is straightforward. 
The numerator t(q;p) comes from the graviton vertices and from the numerator of the free 
propagator. 

Since the integration in ( [AT]) is over all values of q, it is more convenient to make the 
change of variables q «-» — q in all the terms(z <-> —z) so that 



d 3 g 



-ioc+5 



3+5 2n i 



N(z) 



t(q;p) 



Xq • i) % (p • pY 



+ q 



-q 



(A2) 



Factorizing the denominators in ( |A2|) we can write 



, +ioo+5 A 7 

d 3 g / ^N(z) 



z+\q\) (z + k + \p\y 



(A3) 



t(q;p) 



+ q <-> —q 



(z - |g|y (z + k - \p\^ 

The z integration is now readily performed using the Cauchy theorem and closing the contour 
in the right hand side plane where the only poles are located at z — \q\ and z — \p\ — k Q (fc 
is a pure imaginary quantity at this stage of the calculation). In this way we obtain 



d 3 ]im f N(g Q ) t(q; P y 



(A4) 



lim 



9»- 



N(q 



t(g;p) 



U - 1)! go-w-fco dqi- 1 \ (g + h + \ P \y (q ■ q y 
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Performing the change of variables q — > q — k in the second term of ( [^4]) we can write 



I = - d'q 



lim 



i-1 



N(q ) t(q;p) 



i)! ?0 -igi ac 1 vteo + l^ (p-p) j , 



lim 



-i 



A(g ) t{q ,q-k;q + k ,q)* 



(j - 1)! go-kl-fco c^" 1 \{q + k + \q\y (g 2 - \q - k\ 2 ) 



+ q 



(A5) 



Finally, using the property N(qo + &o) — N(qo) and the symmetry g <-> —q we obtain 

1 / N(q ) t(q;p)\ 



d d q 



(i-lYdqt 1 \(qo + \q\y (p ■ p) j J 
1 ( N(q Q ) t(-p;qY 



(A6) 



(j-iy.dqt 1 \(qo + \q\y (p ■ p) 
The special case when i = j = 1 gives the known result 

t(q;p) +t(-p; q) 



q <-> — g 



2? 



(p-p) 



-5 



(A7) 



where the expression inside the bracket is a typical contribution to the on-shell forward 



scattering amplitude. Although the derivatives in the general expression (|A6|) makes it 
much more difficult to be handled, it is straightforward to deal with such kind of expressions 
using a computer algebra program. 



APPENDIX B: 

The Feynman rules are obtained inserting Eq. ( |2.3| ) and the corresponding perturbative 
expansion of the inverse 



-9 9 



-1 



+ 0(k 4 ) 



(Bl) 



into Eq. (|2.1j ). The contributions of order in k yields the graviton propagator given by 
Eq. 
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The third term in Eq. ( |2.1|) yields the following expressions for the ghost propagator and 
the graviton-ghost-ghost vertex 

Vfr\k) = ^, (B2) 



t\j r 

V £°il 2 ^( k ^ fc 2, h) = 2 yl^Ah^hui + fomfovi) - (V^^hu! + V^ 2 ki^)k 2fl3 \ • ( B 3) 



All the graviton self-couplings are generated only from the first term in Eq. (|2.1| ). The 
corresponding Feynman rules for the three- and four-graviton couplings are given respectively 
by the following expressions 

"4 ^2^3 k 3u2 ?7 WM2 Vviv 3 ~ k 2 ■ k 3 ^i^3 Vuius Vij.2V2 + 2k 2 ■ k 3 rj^ lU2 r\ VxV3 ^ 2M3 



4 

+2 k 2 • k 3 Tjfj, lfi2 T]vi/i3 Tlvius 2 k 2fM1 k 3ui ^ 2 a*3 ^2^3 

— /C2 ' ^3 Vvil/2 "I" ^(Ui fel/j V 1*2*2 V^3^3 

+ (symmetrization under (pi Vi), (fi 2 <->• z/ 2 ), (7x3 <-> ^3)) 
+permutations of (/ci, Vi), (k 2 ,[i 2 ,v 2 ), (k 3 ,fi 3 ,u 3 ) , 



(B4) 



H^l/^/^a/M/^ 1 ' ^2, &3, ^4) - 

— —2 fc 3 ^ 2 k^ U2 T}ft lV3 1} VxllA % 4 /i 3 + ^3/i 2 ^ 4 ^2 ^1^3 ^1^3 ^4^4 
k 3 ■ ^4 VfJ-l/J-i VviV2 V/J-4/*4 V/J.2V3 4 ^3^4 ^4^3 V^l^'i Vvi"2 VfJ.2^4 

+2 /C3 • /C4 Tj^ 1U3 Tj ulU2 Tj^ A ^ 3 V/J.2^4 k 3 ■ k^ Tj^^g T\ UXV3 ^ 4( u 2 ^1/2/14 

+2 • Tj^ lfJi3 Tj Vl ^ 4 Vv2^4 ^3^ 2 kiv 2 Vpiw Vvim V^^s ( B ^) 

— /C3 • /C4 ?7^ 1 / J4 ^i^iaj V/J-2H4 Vw3 ~ ^ &3/i 2 ^4y 2 VnxfJ-3 ^1M4 Vv3lM, 
+2 &3 ■ ^4 T]^^ Tj UlU2 ^4^2 ^3/M 

+ (symmetrization under (/ii <-» z/i), (/i 2 <-> z/ 2 ), (^3 <-» 2/3), (^4 <-> ^4)) 
+permutations of (fci, jtii, 1/1), (k 2 ,fi 2 ,v 2 ), (k 3 , /i 3 , z/ 3 ), (fc 4 ,/i 4 , /i 4 ) . 

As usual, we have energy-momentum conservation at the vertices, where all momenta are 
defined to be inwards. 
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APPENDIX C: GRAVITATIONAL T HOOFT IDENTITIES 



The imaginary time formalism at finite temperature follows closely the corresponding 
formalism at T = 0. Consequently, the 't Hooft identities at finite T would be similar to the 
ones at T = 0, were it not for the presence of 1-particle tadpole contributions (such terms 
vanish at T = in the dimensional regularization scheme). However, since the tadpole 
terms are proportional to T 4 , they do not affect the identities involving the sub- leading 
contributions. To derive these, we start from the action 

I = J ^xd A yh, v {x)SZ b aP {x - y)h aP {y) + f d 4 xd*yJ^(x)X^ x (x - y) V \y) + ■■■. (CI) 

Here denotes the sub-leading contributions to the graviton 2-point function and X^ u \ 

represents the tensor generated by a gauge transformation of the graviton field which is 
given to lowest order by Eq. ( |2.6| ) in the momentum space. J^ v is an external source, r/ x 
represents the ghost field and • • • stand for terms which are not relevant for our purpose. 
The 't Hooft identity involving the graviton self-energy function is a consequence of the 
BRST invariance of the action /: 

d A x — —, . J 1 . . = (C2) 

To lowest order, Eq. ( |C2| ) is equivalent to the relation Eq. ( |2.5| ). In general, Eq. ( |C2[) 
implies the generalized 't Hooft identity 

X^ a S^" = 0, (C3) 

which can be written to second order as 

Using Eq. (|2.5|) , we see that ( p^) leads immediately to the 't Hooft identity (|2.9|) . 

In order to derive Eq. ( |3.8|) , we shall need to evaluate the tensor X^J X which appears in 
Eq. (P4|) . This tensor may be represented by the diagram shown in Fig. 4, where the ghost- 
graviton-source vertex is given in the Appendix A of ref. |f[E|. Using the forward scattering 



amplitude method, we obtain the following structure for the T 2 contributions to X 
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X ^ K J dn [ 2r] " x + — WW) (C5) 

Using this form and Eq. ( |2.5| ), it is clear that X^J^ is orthogonal to tJ,ua ^ , so that the 
relation ( |3.8| ) follows at once from the identity 
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TABLES 



rpjJ-lVl _ ^2^1 j^2/H _|_ ^M2Ml ^2^1 

rpjXlVl W2 = ^ ^V2 ^1 + yP2 ^1 ) + ^1 ^2 ^1 + ^1 ) 

T£ 1U1 ^ = u U2 u Ul 

rp/J-lV! p,2V2 _ ^2^2 ^Ml^l 

Tf 1 ^ = U^ 1 U U1 T]^ 2 + U^ 2 U V2 rf 1 " 1 

rp^lVl 112V2 _ U U2 ^M2M1 _|_ k ^l <qH2Vl^ _|_ fcV2 (yV\ ^2^1 _|_ <qH2Vl^ 

+ u^ 2 {k Ul rf 2 ^ + A; Ml rf 2VX ) + A;' 42 (V 1 rf 2 ^ 1 + u^ 1 rf 2VX ) 
Tf 1V1 ^ = k Ul u^ 2 u U2 u^ 1 + k^ 1 u^ 2 u V2 u Ul + k U2 u^ 2 u^ 1 u Ul + k^ 2 u V2 u Ul 
jigivi W2 = k U2 k vi ^2Mi + k U2 k^ 1 7]^ 2V1 + k^ 2 k Ul r] U2f11 + k^ 2 k^ 1 7] V2U1 

T gM V2V2 = ^ k U! yP2 U U2 + ^2 ^2 tfl yVl 

C >* 2U2 = (k U2 u^ 2 + k» 2 u U2 ) (k Ul + k^ u Ul ) 

T Ml^l M2^2 = k V2 k »l k Vl yP2 + £M2 fcjil fc^l U V2 + £M2 ^2 £.^1 ^Ml _|_ £M2 fc« fcMl ^1 

JiVl"! M2^2 _ £/i 2 k V2 k Hl k "\ 

rptg"! V2V2 = ^2 + kW 2 k V2 ^IVI 

rpVlVl V2V2 = ^ ^ + ^2V2 + (£^2 M M2 + £M2 ^2) 

TABLE I. A basis of 14 independent tensors 
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i 








1 

2 

3 

4 
5 

6 

7 

8 

9 

10 

11 
12 
13 
14 


1 

2 

1 k 2 

2 (k-u) 2 -k 2 
1 (A: 2 ) 2 




1 A; 2 

2 (k ■ u) 2 - k 2 





3 (k 2 ) 2 


2 ((k ■ u) 2 - k 2 ) 2 
1 

~2 

1 /fc 2 

2 (k ■ u) 2 - k 2 

1 k ■ u 

2 (k ■ uf - k 2 
1 k 2 k ■ u 


((fe-«) 2 -fc 2 ) 2 

n 
u 



1 • u 

2 (fe • n) 2 - k 2 

k 2 k ■ u 


2 ((k ■ u) 2 - k 2 ) 2 
1 

1 I 2 

2 (k ■ uf - k 2 



3 k 2 k ■ u 


2 ((A;-u) 2 -A; 2 ) 2 
1 2(k-uf-k 2 


2 9 

((fe • u) 2 - A; 2 ) 

1 (k ■ u) 2 

2 fe 2 ((jfe ■ u) 2 - k 2 ) 


2 ((k ■ u) 2 - k 2 f 


1 2(k-uf + k 2 


2 ((fc. u )2 -fe 2 ) 2 
1 k 2 


((A; • u) 2 - A; 2 ) 2 
1 3(k-u) 2 + k 2 


2 ((k-u) 2 -k 2 f 

3 (k ■ uf 


2 ((k-uf-k 2 ) 2 
1 k ■ u 


2 ((fc-u) 2 -fe 2 ) 2 
((fc • u) 2 + k 2 ) k ■ u 


2 ((k ■ u) 2 - k 2 f 
1 (2(k-uf + k 2 )k-u 


2 ((& ■ u) 2 - A; 2 ) 2 
1 ((* ■ u) 2 - k 2 ) 

~\ ((k-uf-k 2 )- 1 

1 k ■ u 

2 (k ■ u) 2 - k 2 


k 2 ((k-u) 2 -k 2 f 
(b i,\ 2 


2 k 2 ((k-u) 2 -k 2 f 

1 A(h . oi\ A .4- A(b . oi\ 2 h 2 -I- (b 2 \ 2 


k 2 {(k-u) 2 -k 2 f 





6 (k 2 ) 2 ((k ■ u) 2 - k 2 f 
1 2(fc-n) 2 + fe 2 

6 k 2 ((k-u) 2 -k 2 ) 

1 fc ■ u 

2 (k ■ u) 2 - k 2 



TABLE II. Components of the transverse traceless tensors 
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FIGURES 




(a) (b) 

FIG. 1. Diagrams contributing to the one-graviton function. Wavy lines represent gravitons 
and dashed lines represent ghosts 




FIG. 2. Diagrams contributing to the two-graviton function. Wavy lines represent gravitons 
and dashed lines represent ghosts 
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(0 i .o 




FIG. 3. The dispersion relations for the modes A, B and C in units of ( 16^(7/?) ^ 2 for real 
frequencies and wave vectors. The dashed lines stand for the leading T 4 contributions and the full 
lines represent the inclusion of the T 2 sub-leading corrections for GT 2 = 0.01. The light-cone is 
represented by the diagonal. 



q + k 



H v 




X 

— 

k 



FIG. 4. The source-ghost Feynman diagram. The full/wave line on the left represents the 
external source. 
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